Quantum Dot and Hole Formation in Sputter Erosion 
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Recently it was experimentally demonstrated that sputtering under normal incidence leads to 
the formation of spatially ordered uniform nanoscale islands or holes. Here we show that these 
nanostructures have inherently nonlinear origin, first appearing when the nonlinear terms start to 
dominate the surface dynamics. Depending on the sign of the nonlinear terms, determined by the 
shape of the collision cascade, the surface can develop regular islands or holes with identical dynam- 
ical features, and while the size of these nanostructures is independent of flux and temperature, it 
can be modified by tuning the ion energy. 
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The fabrication and physical properties of quantum 
dots (QDs) are topics of high current interest due to 
their applications in optical devices and as potential 
construction blocks of novel computer architectures [Q. 
However, despite the high interest in the subject, the 
methods available for the fabrication of such dots are 
rather limited. Lithographic techniques, while under- 
going rapid improvements in resolution still cannot 
produce small and dense enough dots necessary for de- 
vice applications. While much research has focused on 
self-assembled QD formation ||^, that generates dots 
through the unique combination of strain and growth ki- 
netics, these techniques offered relatively uniform islands 
only for a few material combinations and still have to 
find their way into devices. Consequently, there is con- 
tinued high demand for alternative methods that would 
allow low cost and efficient mass fabrication of QDs. In 
the light of these technological driving forces, the recent 
demonstration that sputter erosion can lead to uniform 
nanoscale islands, that exhibit quantum confinement, 
will undoubtly capture the interest of the scientific com- 
munity [Q . Ion beam sputtering has long been a leading 
candidate for surface patterning. While ripple formation 
on sputter eroded surfaces has been observed already in 
the 70s in the last decade much work has been de- 
voted to understand both the experimental and theoreti- 
cal aspects of this fascinating self-organized phenomena. 
However, most experiments have focused on off-normal 
incidence, that, by breaking the symmetry along the 
surface, leads to anisotropic structures, such as ripples. 
While theoretically it was expected that under normal 
incidence the ripples should be replaced by some periodic 
cellular structures, such surface features have not been 
observed experimentally. Recently, two groups have ob- 
tained simultaneous advances in this direction. Facsko et 
al., investigating low-energy normal incident Ar+ sput- 
tering of GaSb (100) surfaces observed that as erosion 
proceeds, nanoscale islands appear on the surface, that 
are remarkably well ordered, and have a uniform size 
distribution. On the other hand, recent experiments of 
Ar^" sputtering of Cu(llO), and Ne+ sputtering of Ag 



(001) under normal incidence lead to relatively uniform 
depressions or holes These experiments, while rep- 
resent significant technological breakthroughs, raise a 
number of questions regarding the mechanisms responsi- 
ble for the formation of these nanostructures. While it is 
tempting to interpret these structures as periodic pertur- 
bations expected by the linear theory of sputtering j7j , a 
careful analysis of the experimental results indicates that 
such approach is less than satisfactory. To mention only 
a few discrepancies, Facsko et al. reported that the ob- 
tained island size is independent of temperature, while, 
according to the linear theory their size should decrease 
exponentially with T Also, the linear theory would 
predict at best a cellular structure displaying a square 
lattice, in contrast with the hexagonal ordering observed 
in experiments 

In this paper we present the first detailed theory ad- 
dressing the formation of sputter- induced QDs and holes. 
We demonstrate that these structures are the result of 
inherently nonlinear phenomena, and thus they cannot 
be accounted for in the context of the linear theory. We 
find that regular nanostructures first appear as the non- 
linear terms become relevant, allowing us to predict the 
characteristic time, t, necessary for their formation, and 
calculate the dependence of r on the physical parameters 
characterizing the ion bombardment process. Further- 
more, we show that the QDs and holes, observed in 
different experiments, are governed by the same physical 
phenomena, the difference between them coming from 
the shape of the ion cascade characterizing the interac- 
tion of the bombarding ions with the substrate. Finally, 
we predict that the size of these islands and holes, while 
independent of flux and temperature, depends on the ion 
energy. 

A particularly successful description of the morpho- 
logical evolution of sputter eroded surfaces has been 
proposed by Bradley and Harper (BH) Q, based on 
Sigmund's theory of sputtering ||^, predicting that the 
height h{x, y, t) of the eroded surface is described by the 



linear equation 
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where Vx and Vy are effective surface tensions, generated 
by tlie erosion process and K is tire surface diffusion con- 
stant induced by thermal diffusion. The balance of the 
unstable erosion term {—\v\d'^h) and the surface diffusion 
term {—Kd^h) acting to smooth the surface, generates 
ripples with wavelength 



= 2n^2K/\u~\, 



(2) 



where i refers to the direction {x or y) along which 
{fx or Vy) is the largest. While successful in predicting 
the ripple wavelength and orientation ||], the linear the- 
ory fails to explain a number of experimental features, 
such as the saturation of the ripple amplitude pO|-p^, 
or the appearance of kinetic roughening HJlsll . To ad- 
dress these questions it has been proposed |16| that these 
shortcomings can be cured by the inclusion of nonlinear 
terms and noise, derived from the Sigmund's theory of 
sputtering. Consequently, Eq. (1) has to be replaced by 
the nonlinear equation, 

dth = y^dlh + vydlh - Dxdth - Dyd^h - Dxyd^d^h 



(3) 



where and Xy describe the tilt-dependent erosion rate, 
depending on flux / and the penetration depth a, and 
ri(x,y,t) is an uncorrelated white noise with zero mean, 
mimicking the randomness resulting from the stochastic 
nature of ion arrival to the surface p^ , [l7t . Furthermore, 
the sputtering process always generates ion induced effec- 
tive diffusion, that together with thermal diffusion, are 
incorporated in D^, Dy, and Dxy WA - Under normal 
incidence, the coefficients in (3) are isotropic, given by 
0@ 



Uy = -faal/2al 



D = D, 
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where = a//i and = a/a and /i and ct, defined in 
Fig. la, characterize the shape of the collision cascade of 
the bombarding ion. 

To investigate the origin of the QDs and the dynamics 
of QD formation under normal incident ion sputter- 
ing, we integrated numerically the continuum equation 
(3), using standard discretization techniques p^ , and 
isotropic coefficients as expected for normal incidence. 
We choose a temporal increment At — 0.01 and impose 
periodic boundary conditions h{x, y, t) — h{x + L, y, t) — 



h{x,y + L,t) where L x L is the size of the substrate. 
To improve the uniformity of the QDs, the numerical 
simulations were carried out without noise (jj — 0), using 
instead a random initial surface configuration. However, 
we repeated the simulations for the noisy case as well, 
finding that the size uniformity of QDs is enhanced as 
the noise amplitude decreases. 

Our main result, presenting the morphology of the 
ion sputtered surface at three different stages of their 
time evolution, is shown in Fig. 2. Let us first concen- 
trate on the A > case (upper panels in Fig. 2). In 
the early stages of the sputtering process the surface is 
dominated by small, wavy perturbations (Fig. 2a) gener- 
ated by the interplay between the ion induced instability 
and surface relaxation. However, since the system is 
isotropic in the {x, y) plane, these ripple precursors are 
oriented randomly, generating short wormlike morpholo- 
gies on the surface. After some characteristic time, r, 
these structures turn into isolated but closely packed 
islands, reminiscent of the QDs reported experimentally 
(Fig. 2b). Note that upon a closer inspection one can 
observe the emergence of hexagonal order in the island 
positions. As the sputtering proceeds, while the islands 
do not disappear, the supporting surface develops a rough 
profile, destroying the overall uniformity of the islands 
(Fig. 2c). A similar scenario is observed for A < 0, the 
only difference being that now the islands are replaced 
by holes (Figs. 2d-f), reminiscent of the morphologies 
observed experimentally by Rusponi et al. The first 
conclusion we can draw from these results is that the 
development of QDs and holes is governed by the same 
underlying physical phenomena, the only difference being 
that for QDs we have A > 0, and for holes A < 0. Indeed, 
this morphological change is expected from the nonlin- 
ear continuum theory, Eq. (3) being symmetric under 
the simultaneous transformation A — > —A and h — > —h, 
indicating that changing the sign of A does not affect the 
dynamics of the surface evolution, but simply turns the 
islands into mirrored holes. Since, according to (6) the 
sign of A is determined only by the relative magnitude 
of flcr and a^, whether islands or holes appear is deter- 
mined by the shape of the collision cascade, shown in 
Fig. la. Consequently, using (6) we can draw a phase 
diagram in terms of the reduced penetration depths a„ 
and that separates the regions displaying QDs versus 
holes (Fig. lb). These results also indicate that the QDs 
and holes are inherently nonlinear objects, that can be 
explained only by the nonlinear theory, since, should 
the linear terms be responsible for their formation, the 
surface morphology should not depend on the sign of 
A (indeed, Eq. (1) has a full h —h symmetry). In 
the following, we investigate the dynamics of QD and 
hole formation, providing further proof of their common 
nonlinear origin. 

The crossover behavior from the linear to the nonlin- 
ear regimes can be monitored through the surface width. 
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W'^{L,t) = Y^J2x,yh'^i^:y^t)-^'^- As Fig. 3 shows, this 
quantity exhibits a sharp transition at a characteristic 
time t: for t < t, the width W increases exponentially 
as ~ exp{i't / i'^) , while for t > r, still increases 
but at a considerably smaller rate than an exponential 
po| . The crossover time is given by r ^ cr^//a in terms 
of the experimental parameters pO| , indicating its de- 
pendence on the flux and ion beam energy. Correlating 
these results with the observed surface morphologies, we 
find that the QDs first appear at t « r. Indeed, in Fig. 3 
we marked with arrows the time when the morphologies 
in Fig. 2a-c were recorded, indicating that no QDs exist 
before r (Fig. 2a), but they are fully developed at r 
(Fig. 2b), and their uniformity rapidly vanishes after r 
(Fig. 2c). 



coincides with £. Indeed, we determined the diameter 
of the QDs in simulations for i/ = —0.6169, D = 2 and 
A = 1, observing approximately 16 QDs aligned along 
the side of the system of the size 256 x 256, in agreement 
with the prediction L/i « 16. Using (4) and (5), we find 
that I = v^TT/i, i.e., from the average separation of the 
islands one can determine the size of the horizontal 
width of the collision cascade (Fig. la). Furthermore, 
since typically we have ~ a ~ e^™, when e is the ion 
energy and m is an constant that weakly depend on e 
(m = 1/2 for e w 10 ~ 100 keV) [|, we predict that 
one can tune the size of the QDs by changing the ion 
energy e, while the size is independent of the flux and 
the temperature. 



Based on the results presented in Figs. 2 and 3, the 
following scenario emerges for QD and hole formation. 
In the early stages of the erosion process the linear the- 
ory correctly describes the surface evolution, and thus 
we observe the cell structure predicted by the BH theory. 
However, as time increases, the nonlinear terms turn on 
breaking the up-down symmetry of the surface. Such 
delayed effect of the nonlinear terms is a well known fea- 
ture of the class of surface evolution equations which (3) 
belongs to ||2l[] . The sign of A determines whether QDs or 
holes form, these structures appearing at a characteristic 
time T, that depends on the flux and the ion energy, and 
their characteristic size is given by the ripple period ^ 
(Eq. (2)). In turn, the sign of A is determined only by the 
shape of the collision cascade (Fig. 1). As time increases 
beyond r, the nonlinear terms lead to kinetic roughening 
of the surface at large length scales [^21|, and while 
the QDs and holes do not disappear, the substrate on 
which they exist becomes rough, destroying the overall 
island/hole uniformity and ordering. Indeed, as Fig. 3b 
shows, the island height distribution has a larger width 
after r than right at r. Thus, ordered and uniform QDs 
and holes can be obtained only in the early stages of 
the nonlinear regime, more precisely, at the crossover 
between the linear and nonlinear regimes emerging at 
the characteristic time r. 

Experimentally the crossover time r can be estimated 
through monitoring the evolution of the ripple amplitude 
using in situ spectroscopic techniques pl| . Our results 
indicate that in order to obtain uniform islands one needs 
to stop the sputtering process just when the amplitude 
of the eroded surface saturates: before saturation the 
QDs have not developed yet, while after saturation their 
uniformity disappears. Since the mechanism generating 
the QDs is the ion induced instability, balanced by ion 
induced diffusion, their typical separation is given by the 
most unstable wavelength, I = 2TTyj2D /v. The simula- 
tions indicate that the islands (or holes) first appear as 
small surface perturbations. However, as their amplitude 
increases, around r they get in contact with each other, 
and their growth saturates. Thus their final diameter 



Finally, the conclusion i ^ fj, holds only when ion 
induced effective diffusion JlTj is the main relaxation 
mechanism. While the temperature independent QD 
size reported by Facsko et al. [|j indicates that for GaSb 
this is the main relaxation mechanism, at high temper- 
atures or in othe r sy stems thermal diffusion could be 
more relevant [|l^,|l8[^. Note, however, that the sce- 
nario presented above for QD and hole formation, is 
not conditional on ion induced diffusion. Should ther- 
mal diffusion be the dominating relaxation mechanism 
(certainly true in any system at high temperatures) it 
will change only our prediction for £ and r. In this case 
we expect £ ^ 2-k^2I:)q exp(-£;/2fcr)/i/, thus the QD 
size will depend exponentially on T, and we also have 
i ^ e^""^/^, and I ^ f^^^^- For the crossover time we 
obtain r ~ exp{-2E / kT) / , r - e-\ and r ~ f~^. 
However, the phase diagram (Fig. lb) and the expected 
dynamical evolution (Figs. 2 and 3) will not be sensitive 
to the nature of the relaxation mechanism. 



In conclusion, we demonstrated that the nonlinear the- 
ory can fully explain the recently observed quantum dots 
and holes generated by sputter erosion. We showed that 
these nanostructures appear only when the nonlinear 
terms become effective, stabilizing the surface amplitude, 
and that the difference in the islands and holes comes 
in the sign of the nonlinear terms. Furthermore, we are 
able to predict the dependence of the characteristic time 
for QD formation, and find that their saturation size is 
given by the size of the collision cascade. We believe that 
these results will lead to a better understanding of the 
formation and evolution of these fascinating nanostruc- 
tures, thus they will guide further experiments aiming to 
better control these systems. 
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FIG. 1. (a) The characteristics of the collision cascade gen- 
erated by an ion. According to Sigmund's theory of sputter- 
ing [8] , an ion penetrates the surface until a peneration depth 
a, where it radiates its energy out, such that the kinetic en- 
ergy generated damage, decays exponentially with the dis- 
tance from O. For most systems the energy distribution is 
anisotropic, being characterized by the characteristic decay 
lengthscales a and fj, along or perpendicular to the ion direc- 
tion, respectively, (b) Phase diagram showing the parameter 
regimes corresponding to island (QD, A > 0) and hole (A < 0) 
formation. 
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FIG. 2. (a)-(c) Surface morphologies predicted by Eq. (3), 
for A = 1 at different stages of surface evolution. The pic- 
tures correspond to (a) t =4.0, (b) 5.8, (c) 8.0x10*. (d)-(f) 
The same as in (a)-(c), but for A = —1. In all cases we used 
V = 0.6169, D — 2, and system size 256 x 256. 
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FIG. 3. (a) Time evolution of the surface width for the 
parameters v = —0.6169, D = 2, A = 1 or —1, and system size 
256 X 256. The crossover time is estimated to be r ~ 5.8 x 10'* . 
The arrows denote the moments when the surface snapshots 
shown in Fig. 2 were recorded, (b) Island height distribution 
right at the crossover time t = r (black columns), and at 
a later time, t=8.0xl0* (gray columns), indicating that the 
width of the distribution is smaller at the crossover time, and 
broadens as we go beyond r. 
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